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We derive the equations that describe adsorption of diffusing particles onto a surface followed by 
additional surface kinetic steps before being transported across the interface. Multistage surface 
kinetics occurs during membrane protein insertion, cell signaling, and the infection of cells by virus 
particles. For example, viral entry into healthy cells is possible only after a series of receptor 
and coreceptor binding events occur at the cellular surface. We couple the diffusion of particles 
in the bulk phase with the multistage surface kinetics and derive an effective, integro-differential 
boundary condition that contains a memory kernel embodying the delay induced by the surface 
reactions. This boundary condition takes the form of a singular perturbation problem in the limit 
where particle-surface interactions are short-ranged. Moreover, depending on the surface kinetics, 
the delay kernel induces a nonmonotonic, transient replenishment of the bulk particle concentration 
near the interface. The approach generalizes that of Ward and Tordai [1] and Diamant and Andelman 
[2] to include surface kinetics, giving rise to qualitatively new behaviors. Our analysis suggests a 
simple scheme by which stochastic surface reactions may be coupled to deterministic bulk diffusion. 

PACS numbers: 68.43. +h, 87.68. +z, 68.47.Pe, 68.03.+g 



I. INTRODUCTION 

The kinetics of surface particle adsorption and of 
transport through interfaces play a key role in sur- 
factant phenomena [3, 4], membrane biology and cell 
signaling [5-10], marine layer oceanography [11], and 
other biological and chemical processes. Particle ad- 
sorption may fundamentally alter the physical and 
chemical properties of the interface, and it is crucial 
to understand both equilibrium and dynamical prop- 
erties of the adsorbed layers [1, 3, 4, 6]. In the seminal 
work of Ward and Tordai [1], a bulk phase acting as a 
reservoir of particles is physically limited by an empty 
surface onto which the particles can adsorb. Particles 
are assumed to lower their free energy with respect 
to the bulk phase by irreversibly and instantaneously 
adsorbing onto the interface. Under these conditions, 
the total concentration of adsorbed particles may be 
estimated in relation to measurable interfacial prop- 
erties, such as the dynamic surface tension. Several 
applications, extensions and alternate approaches to 
this work have been proposed [2, 12]. In particular, 
adsorption dynamics in the Ward- Tordai setting can 
be rederived through a free energy approach [2], al- 
lowing for the inclusion of ionic surfactant effects and 
electrostatic interactions. 

In many biochemical systems, the complete adsorp- 
tion of a particle arriving from the bulk requires a 
series of auxiliary transformations at the surface be- 
fore the particle can be successfully incorporated, or 
'fused' into the surface. These intermediate steps 
gives rise to a lag-time in the complete adsorption 
process. For example, the incorporation of emulsify- 



ing proteins onto an air-water interface may be de- 
layed by the unfolding of the polypeptide at the inter- 
face [7]. Adsorption of proteins on polymer-grafted 
interfaces, such as the glycocalyx layer of vascular 
endothelial cells, is also delayed due to the progres- 
sive insertion of the protein through the polymer 
brush[13, 14] . Kinetic delays have also been observed 
in the adsorption of the hemagglutinin glycoprotein 
(HA) of the influenza virus as it enters target host 
cellular membranes [15]. The mechanisms underlying 
this delay are not known in detail but are believed to 
involve conformational changes of HA molecules into 
fusion enabling complexes, mediated by the presence 
of binding receptors and coreceptors on the target 
cell membrane [15-17]. Similarly, the incorporation 
of an HIV particle into a T-cell or a macrophage is 
possible only after the gpl20 glycoprotein of the HIV 
virus membrane recognizes and binds to the target 
cell surface receptor CD4, and subsequently to other 
coreceptors such as CCR5 or CXCR4. As in the case 
of HA and influenza, the exact number of gpl20- 
bound receptors and coreceptors required for HIV 
particle fusion is yet unknown and might depend on 
gpl20 conformations and receptor/coreceptor bind- 
ing cooperativity [18, 19]. The complex nature of 
surface biochemistry makes quantitative kinetic mea- 
surements challenging. Recently, the binding kinetics 
of the CD4 cellular receptor to the gpl20 HIV lig- 
and have been measured under different experimental 
conditions yielding widely different dissociation rates 
[20, 21]. In this work, we will provide a quantitative 
framework that can be used to better understand the 
experimentally observed lag-times in surface kinetics 
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FIG. 1: A schematic of a typical adsorption experiment 
on a confluent monolayer of target cells. The location of 
the interface is labelled by S. Particles such as viruses 
are spread over the cell layer in a thin supernatant film. 
Inset: After initial nonspecific viral adsorption on the 
supernatant-cell interface, cellular receptors and corecep- 
tors bind to the virus via a certain stoichiometry, forming 
fusion intermediates IY The subsurface layer of thickness 
I, the subsurface concentration n s (t), and the adsorbed 
species Ti(t) are discussed in the text. 



chemical systems [23]. In Sect. 4 we particularize our 
surface reaction scheme to a specific Markov process 
chain and evaluate all physically relevant quantities. 

A. Bulk Diffusion 

In the continuum limit, the density of parti- 
cles n(r,t) in the bulk phase obeys the convection- 
diffusion equation 

d 1 

= V ■ [DVn] + — V • [D nVU] , (1) 

where D(r) and U (r) are the local diffusion coefficient 
and potential of mean force, respectively, and fc^T is 
the thermal energy. Spatial variation of D(r) and 
U(r) may arise from interactions with the interface 
as shown in Fig. 3. Boundary conditions are typically 
applied at the mathematical surface onto which the 
particles adsorb or reflect. By balancing the diffusive 
flux just above this mathematical interface with the 
particle rate of insertion into the interface, a mixed 
boundary condition arises 



phenomena that involve multistage surface chemistry. 

In particular, we will explicitly consider intermedi- 
ate, reversible steps for surface binding in the Ward- 
Tordai formalism, deriving an effective boundary 
condition to complement the bulk diffusion process. 
Chemical transitions among the surface species will 
introduce memory terms in the boundary conditions 
for the bulk concentration. Our analysis can be read- 
ily applied to the titration of replication-incompetent 
virus via a colony formation assay [22] as shown in 
Fig.l. 



II. MODEL EQUATIONS 

In this section, we motivate and derive the equa- 
tions coupling bulk diffusion to surface layer evolu- 
tion. We consider a general, linear reaction scheme 
to describe the multistcp surface reaction dynamics. 
Effective boundary conditions for diffusion from the 
bulk are derived in Sect. 3. As we shall discuss in 
detail, we are able to embody the response of the 
adsorbing particle system to the existence of inter- 
mediate chemical steps at the surface, into a unique 
delay kernel regulating the boundary dynamics. All 
microscopic details stemming from the surface dy- 
namics, no matter how complicated, are contained in 
the derived memory kernel. Our approach includes 
ligand rcbinding to surfaces, found to be important 
for analyzing surface plasmon resonance assays of bio- 



D(r)h- Vn(r,i) =yn(r,t), reS. (2) 

Here, S denotes the substrate; its normal direction is 
n. The parameter 7, which has the physical units of 
speed, is proportional the probability / (often called 
the accommodation coefficient [24] or sticking proba- 
bility [25, 26]) that a particle is adsorbed into the 
mathematical interface upon collision. We define 
7 = 7o/ such that in the limit 70 — > 00 and / 7^ 0, 
Eq. 2 is equivalent to n(r G S) = 0, an absorbing 
boundary condition. A reflecting boundary condi- 
tion, n • Vn = 0, arises when / = 0. Equations 1 
and 2 are commonly used to model simple diffusion- 
adsorption processes at surfaces. 

B. Surface Reactions 

In many applications, particles at an interface un- 
dergo chemical or physical modifications that control 
for example, surface reactivity, surface tension [1-4], 
and conductivity [25]. Biological examples include 
tissue factor initiated coagulation reactions and vi- 
ral entry. Coagulation factors must work their way 
through the glycocalyx layer before they can be en- 
zymatically primed by the membrane-bound tissue 
factors [27]. Entry of viruses, such as HIV, into 
cells require the binding of membrane-bound recep- 
tors and coreceptors before fusion with the target cell 
can occur. All of these processes can be thought of 
as reactions at the membrane surface. Immediately 
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FIG. 2: Surface reaction scheme among the intermediate 
states rV The source species Fi is supplied by the bulk 
surface concentration at the interface, n s . In the context 
of virus recognition and infection, the intermediate steps 
label various numbers of receptors or coreceptors associ- 
ated with the surface bound virus particle. For example, 
Yi may denote the surface concentration of HIV particles 
with i receptors and coreceptors attached. This catenary 
model could also represent successive degrees of insertion 
of an absorbing species through a polymer brush or gly- 
cocalyx coated interface. 



after adsorption from the bulk, the surface particle 
concentration, whether of coagulation factors or of 
virus particles, is denoted by T\. For example, in 
the case of viruses, we can identify the Ti state as 
being that of a virus bound to i = 1 CD4 surface 
receptor. The initially adsorbed species can then ki- 
nctically evolve into the other species Ti representing 
virus particles with i > 1 bound receptors or corecep- 
tors. The kinetics among the N surface species fol- 
lows the linear rate equation d t T = Mr+F(i), where 
r = (ri, T2, ■ ■ ■ , Tjv), M is the transition matrix 
among the N surface states, and F = (F, 0,...,0) 
is the source of the first, originating source species 
Ti coming from the bulk. 

Figure 2 illustrates a simple example of a linear 
surface reaction scheme that can be described by the 
above linear rate equation. In this case the reaction 
matrix M is tridiagonal. General reaction matrices 
can also be analyzed since our results depend only on 
the eigenvalues and eigenvectors of M. 
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FIG. 3: Interaction potential between diffusing species 
and the interface. The subsurface layer is defined by ei- 
ther the range of the interaction potential lint, or the 
mean free path £ m f P , depending upon which is larger. 

must solve Eq. 1 with a nonuniform U(r), and possi- 
bly a nonuniform D(r), to within a distance £ = £ m f p 
of the interface. For the choice £ = £ m f P , the adsorp- 
tion velocity 70 can be approximated by the thermal 
velocity vt such that 7 = 70/ ~ vxf- The value of 
the bulk density at the boundary S + £h is defined 
as the sublayer density: n(r = S + £h,t) = n s (t). 
The equation for the rate of change of the number 
per area of molecules in the thin layer, d(£n s (t))/dt 
can be obtained by balancing the latter with the dif- 
fusive flux into the layer D(r) n • Vn(S + M, t), the 
adsorption into the surface concentration Ti, and the 
spontaneous desorption from the initially adsorbed 
species Ti occurring at rate qi. The complete set of 
equations coupling the variables n(r, t), n s (t), and Ti 
is thus 



C. Surface Layer 

Because the surface densities Ti carry units of num- 
ber per area, and the bulk densities are expressed by 
number per volume, any kinetic parameter linking 
bulk source concentrations to those at the interface 
must introduce a physical length scale. Diamant and 
Andclman [2] have introduced the sublayer thickness 
as a mathematical step coupling the bulk density to 
surface density. Here, we physically motivate this 
"surface layer" and the associated transport. Let us 
thus introduce a thin layer of thickness £ near the 
surface, in which the particle density is denoted n s (t) 
and is still expressed in units of number per volume. 

The continuum approximation Eq. 1 breaks down 
when resolving the transport within distances of a few 
mean free paths. If we identify the sublayer thickness 
£ with the mean free path £ m f P , as shown in Fig. 3, we 



On 
~di 

' dt 
dT 
dt 



= V-[UVn] + -— 
= -F + £>n-Vn 



V-[DnVU], (3) 
qi fi, (4) 



Mr + F, 



(F,0,0,...,0). (5) 



Here, n\ r= s+i n = n Sl and F = F[n Sl {Ti}] is the flux 
of the surface concentration n s into the incipiently 
adsorbed species T\. This functional may depend on 
interactions among the adsorbed species Ti , including 
cooperative or crowding effects, and may be modeled 
using free energies and chemical potential differences 
between the bulk and surface [2, 13]. 

A further simplification can be introduced by defin- 
ing a different sublayer thickness £ = £i n t , where £i n t 
is the typical range of the particle-surface interaction 
as shown in Fig. 3. In this case, at least at a distance 
£i nt from the interface, D(r) is constant, U(r) is zero 
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and Eq. 1 is approximated by the standard diffusion 
equation 



dn(r, t) 
dt 



D\7 2 n(r,t). 



(6) 



All effects of the potential of mean force U(r) and 
spatially varying D(r) are now subsumed into an ef- 
fective source F[n s , {1^}]. This is consistent with all 
previous treatments [1, 2, 4] in which transport in 
the bulk phase was described by simple diffusion with 
uniform D and U. Provided F[n s ,ri] is independent 
of I\, equations 1-5 can be explicitly solved in simple 
geometries. For low surface densities I\ such that 
additional adsorption is not hindered by steric ex- 
clusion, one can assume a form F[n s ] independent 
of surface concentrations IV For £ = tint, 7 is now 
interpreted as an effective adsorption coefficient al- 
lowing Eq. 1 to be replaced by Eq. 6, and simplifying 
the bulk concentration equation. 



III. ANALYSIS 

Following the original work of Ward and Tordai, 
subsequent studies on adsorption and dynamic sur- 
face tension measurements [1, 4] eliminate the bulk 
density at the interface in Eq. 4 to yield two coupled 
integro-differential equations for n s and T\ which 
must be numerically sclf-consistently solved. Here, 
we solve the linear Eqs. 5 independently from the bulk 
densities, but with a source term F[n s ] that connects 
the surface concentrations T,; with the bulk concen- 
tration n(r,t). If F[n s ] is independent of I\, the 
explicit solution to Eqs. 5 can be found by evaluat- 
ing the eigenvalues Xj and corresponding eigenvectors 
v- 7 of the chemical transition matrix M. Denoting 
the similarity transform matrix Vkj = v J k such that 
VMV -1 =diag(Aj), the surface densities are 



A' 



m) = E 



x.i 



N 

E 

3=1 



(0)e 



e^^") F[n s {t')]dt' , 



(7) 

where T(0) are the intermediate surface concentra- 
tions at t — 0. If there are no spontaneous sources of 
the surface cell intermediates, all eigenvalues Xj < 0. 
From Eq. 7, in the case T(0) = 0, T(t) is proportional 
to F[n s (t)). Upon substituting Ti(t) 

By setting k = 1 in Eq. 7, we substitute Ti(t) into 
Eq. 4, and find a concise description of the diffusion- 
adsorption process: 



dn(r,t) 
dt 



' dt 



DV 2 n(r,t), n{S+£n,t) =n s (t), (8) 



Dn-Vn(r,t) 



JV 



K(t-t')F[n s (t')]dt', 
(9) 



K(t)=S(t)-q 1 J2Vr j 1 V jl e^ t , 



(10) 



where K(t) is the kernel constructed from the eigen- 
values and eigenvectors of M. It is composed of an 
instantaneous response - the immediate depletion of 
n s due to adsorption into the Ti surface species - and 
delay terms arising from the surface kinetics of Eq. 5. 
The complete set of equations 8-10 is one of our main 
findings. This result explicitly shows how multistage 
adsorption is modeled by a bulk diffusion equation 
with an nonlinear integro-differential boundary con- 
dition that incorporates the delay arising from the 
multistep kinetics. Under this scheme, all effects of 
surface reactions are incorporated in the kernel K(t). 

Our analysis can be carried further by specifying a 
linear form for the r^-independent source term 



F[n s (i)]= 7 n s (t), 



(11) 



which simply takes the source for the surface concen- 
tration Ti of the first species to be proportional to 
the subsurface concentration. The surfaces densities 
Ti(t) can be found by substituting n s (t), derived from 
Eq. 8-10, into the expression for F[n s (t)} in Eq. 7. 
Note that the boundary condition Eq. 9 contains a 
singular perturbation, and that for times t £/"/, 
the "outer solution" approximation £(dn s /dt) ps 
yields the standard mixed boundary condition Eq. 2 
with an additional memory kernel. Moreover, in the 
linear approximation of Eq. 11, the convolution of 
the delay term in the effective boundary condition 9 
is amenable to analysis by Laplace transforms. 

For simplicity, we will assume a simple one- 
dimensional problem where all quantities vary spa- 
tially only in the direction normal to an infinite, flat 
interface at z = 0. We nondimensionalize all quanti- 
ties by using I as the unit of length, and q^ 1 as the 
unit of time. Henceforth, in all equations, we make 
the replacements z — > z/t, t — > qii, n — ► £ 3 n, n s — ► 
£ s n s , -» £ 2 G i} D £~ 2 D/ qi , and 7 -> 
To render the notation less cumbersome we omit the 
bars from the redefined quantities. In the discus- 
sion that follows, the z, t, n, n s , T^, D, 7 parameters 
are intended as nondimcnsional. Upon taking the 
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Laplace transform in time of the dimcnsionless forms 
of Eqs. 8, 9, and 10, we obtain 



sh(z, s) — tiq = Dd^n(z, s), 



(12) 



sh s (s) — no = Dd z n(z, s) -jK(s)n s (s), (13) 

z— 1 

where 



k( s ) 



(14) 



and no is the initial, dimcnsionless constant bulk and 
sublayer concentration. The general solution to the 
bulk density h(z, s) from Eq. 12 is 



h(z,s) 1 7K(s)exp[-(z-l)i/s/D 



n 



s(s + VsD + yK(s)) 



(15) 



Once the bulk density is derived, all other quantities 
can be found by inverse Laplace transforming n(z, s). 
In the absence of spontaneous sources of the surface 
intermediates, Xj < 0. In this case, it is possible to 
show that n(s, z) only has a simple pole at s = and 
a branch cut on s = (-co, 0]. Performing the integral 
along the latter, we find the exact results 



n(z, t) 



n Q I L(z,u)e- Ut du, (16) 
'o 



and 



N 



r fc (t) = -mY, v v v » 

where 



u + Xj 



■L(l, u)du, 
(17) 



L(z,u) 



1 (u- 1 k{-u))smJ%{z-l) 



u(u - jK(-u)) 2 + Du 1 



1) - 

-LK{-u). 



(18) 



r 



Equations 16-18 are used to numerically compute all 
of our results in the next Section. For complete- 
ness, analytic expressions for asymptotically short 
and long time limits arc derived in the Appendix. 



IV. RESULTS 



We now specify a surface reaction scheme and con- 
struct its delay kernel by using its associated eigenval- 
ues and eigenvectors. For applications such as multi- 
ple receptor binding of the adsorbed species, we con- 
sider a reversible sequential Markov process among N 
chemical intermediates Ti, i = 1, . . . , N. As shown 
in Fig. 2, formation of state IV|_i from state 1^ oc- 
curs at rate pi , while the reverse step occurs at rate 
<7i+i. The final state Tjv is irreversibly annihilated, 
by transport across the membrane, or by irreversible 
reaction, with rate p*. We can then explicitly write 
the dimensional form of Eq. 5, where M is a tridiag- 
onal transition matrix, as 



d£i 
dt 

dT 
~dt 

dJV 
dt 



F[n s , {rj] - [pi + q 1 )T 1 + q 2 T 2l 



-tt = Pt-i r i-i - (li + Pi) r i + 9;+i r *+i 2 < i < N — I 



= -(p* + q N )T N +p N -iT N -i. 



(19) 

In the simplest case where there is only one surface 
intermediate before transport across the interface, 

= 1 and the dimensionless (qi = 1) delay kernel is 
simply K(s) = — (s + p*)/(s + p* + 1). The sublayer 
concentration n s (t) derived from Eq. 16, and the sur- 
face concentration evaluated from Eq. 17, are shown 
in Figs. 4 for various values of 7. 

Let us estimate typical parameter values for viral 
fusion or molecular binding processes. Typical dif- 
fusion constants for viruses of diameter lOOnm and 
in aqueous environments, are D ~ 10 _8 cm 2 /s. Us- 
ing the typical screened electrostatic interaction po- 
tential, t w 10~ 7 cm, we estimate the dimcnsionless 
diffusion coefficient D ~ 10 6 s _1 /(7i. On the other 
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FIG. 4: Surface densities n s (t) for N = 1. (a) The sublayer density n s (t) as a function of time, for various values of 
the dimensionless adsorption rate 7. The other parameters are fixed at D = 10 4 and p* = 1. (b) The sublayer density 
n s (i) for D = 10 4 , and p* = 0.1. Note the bump in concentration imparted by the slow annihilation rate p* . (c) The 
deviation of n s (t) (for 7 = 1000) if the transient term ldn s /dt is neglected. The solid curves are found from the full 
equation 13, while the dashed curves are solutions when the left-hand-side of Eq. 13 is neglected. The deviation occurs 
only at very short times, is independent of p* , and is greatest for smaller diffusion coefficient D. 



hand, typical diameters of small ligand molecules are 
of the order of lnm yielding a nondimcnsional diffu- 
sion constant D ~ 10 8 s~ 1 /qi. The nondimcnsional 
7 estimated using the thermal velocity vt is now 
7 w fvT/(£qi). For virus particles 7 ~ 10 8 s _1 //q r i, 
while for molecular ligands 7 ~ 10 10 s _1 //gi. The 
dissociation rate q\ is highly variable and typically 
falls in the broad range qi ~ 10~ 4 s" 1 to 10 4 s _1 . For 
the gpl20-CD4 interaction, the dissociation has been 
estimated in model systems [20] to be q± ~ 10~ 3 s~ 1 , 
while the detachment rate for mutant viral species 
[21] can be as high as q\ ~ 10 3 s _1 . Lower dissociation 
rates are possible in tighter binding ligand receptor 
pairs such as EGF-receptor [28] where q\ ~ 10 _4 s . 
For other pairs such as P-selectin and its receptors 
[29, 30], qi ~ 0.1s" 1 to 1 s^ 1 . The sticking proba- 
bility / is proportional to the binding probability of 
upon ligand- receptor contact, multiplied by the re- 
ceptor area fraction at the interface. The factor / 
depends on the receptor density, but is typically of 
the order / - 10~ 4 - 10~ 2 . 

In Fig. 4(a) we plot the sublayer density n s as 
a function of time. For p* = 1, Fig. 4(a) shows 
that the sublayer density n s (t) starts at its initial 
value no and decreases with a nondimcnsional rate 
proportional to 7, eventually monotonically reaching 
n s (t — > 00) — > 0. If the annihilation rate p* is de- 
creased, n s may no longer be monotonic. The ob- 
served increase in the surface concentration is due to 
the slow consumption of material at the interface, al- 
lowing some of the material to desorb after being de- 
layed at the interface, rather than irreversibly reach- 
ing the final annihilated or fused state. For example, 



when p* = 0.1, the surface concentration n s (t) first 
decreases but recovers slightly at longer times, before 
ultimately decaying to zero as shown in Fig. 4(b). 

Note that the curves for n s (t) exhibit a transient 
at short times determined by 1/7. Beyond this 
transient, the full solution we plot in Fig. 4(b) re- 
duces to the outer solution, corresponding to setting 
dn s /dt — on the left hand side of Eq. 9. Fig. 4(c) 
explicitly shows different behaviors of the full and 
outer solutions during the transient time. The effect 
of the idn s /dt term is to slow down the initial de- 
crease of n s , particularly for short times within the 
transient defined by 1 /j. The effects of neglecting the 
boundary layer are less pronounced for larger bulk 
diffusivitics D. 

The temporal evolution of n s (t) is strongly depen- 
dent on Ti(t). In fact, the nonmonotonicity of n s (t) 
for small p* shown in Fig. 4(b) arises from the build- 
up of Ti indicated in Fig. 5(a) which can get released 
back into the subsurface layer. For smaller p*, Ti 
reaches larger values. As long as p* > 0, both n s and 
Ti vanish at sufficiently long times. Complete par- 
ticle depletion near the surface occurs in dimensions 
less than two because there is no bounded steady- 
state solution to the diffusion equation and the de- 
pletion zone moves away from the interface for all 
times as shown in Fig. 5(b). Despite free diffusion, 
the bulk is unable to sustain a particle source near 
the surface as is known from classic diffusion theory 
[31]. The replenishment at small annihilation rates p* 
also manifests itself in the bulk. In the case shown in 
Fig. 5(b), as time increases from t = 0.75 to t = 6, the 
bulk concentration near the interface recovers before 
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FIG. 5: (a) The surface concentration for various 

p* with D = 10 4 and 7 = 200. (b) The bulk density 
profiles n(z, t) as a function of position z at times t = 
0.01,0.1,0.75, and 6. 



ultimately decreasing according to Eq. 21. 

For general AT, the eigenvectors and eigenvalues 
must be explicitly computed. In Fig. 6(a), we plot 
n s as a function of N for uniform q = 1 and uni- 
form p. For small p, the surface kinetics is a highly 
biased random walk away from Tn toward Ti, result- 
ing in a larger n s . Both small p and large N, hinder 
the annihilation process and impart a more reflective 
character to the interface. After initial transients, 
both n s and the surface concentrations T.j maintain a 
high level for a long time before dissipating. Larger 
N also effectively trap surface material in the surface 
reservoir IV The relative amounts of Tj for N = 4 
are shown in Fig. 6(b). For the small p = 0.1 used, 
most of the surface density lies in the initial species 
Ti, decreasing in the latter species. 



V. SUMMARY AND CONCLUSIONS 

We derived an effective, integro-differential equa- 
tion for the boundary condition of a simple diffusion 
process. The approach presented differs from the typ- 




60 timer 120 



FIG. 6: Dependence of surface quantities on number of 
steps in the surface reaction scheme, (a) The sublayer 
density n s as N is increased. The initial rapid fall from 
rig/no = 1 is imperceptible on this scale, (b) The surface 
concentrations I\j as a function of time for iV = 4. 



ical Ward and Tordai treatments since we use a linear 
model for the time rate of change of the initially ad- 
sorbed species Ti, rather than eliminating the bulk 
diffusion equation. The effects of intermediate chem- 
ical steps at the boundary are described by a delay 
kernel that can be decomposed using Laplace trans- 
forms. This kernel is an explicit function of the eigen- 
values and eigenvectors of the surface reaction tran- 
sition matrix. Our results suggest that measurement 
of a few quantities, such as fluorescence monitoring of 
the sublayer density [32], can be used to reconstruct 
the principal components of K(t). This approach can 
be used to probe qualitative features of the surface 
kinetics important in modeling cell membrane signal- 
ing and viral infection, where a sequence of chemical 
steps at the surface are required before initiation of 
signaling or viral fusion. In HIV infection, the initial 
adsorption rate would be proportional to the surface 
CD4 concentration, and the subsequent rates in the 
reaction scheme in Fig. 2 would depend on the core- 
ceptor concentrations, their surface mobilities, as well 
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as the effects of cooperative binding [18]. All of these 
physical attributes are encoded in the distribution of 
eigenvalues and eigenvectors of M. 

For simple linear reaction schemes on a flat sur- 
face, we find explicit dependences of the surface con- 
centrations Ti and sublayer concentrations on the 
eigenvalues and eigenvectors of the transition matrix. 
For smaller annihilation rates p*, and at least one 
(N > 1) surface intermediate, we find that the sur- 
face concentration persists and can replenish the bulk 
concentrations n s after its initial decay. The deple- 
tion zone in the bulk can also recover. Delays that 
induce instabilities in dynamical systems have been 
well established [33]. Here, although the delay oc- 
curs in a boundary condition, we observe nonmono- 
tonic behavior arising in the bulk concentrations as 
well. This rebounding effect is also apparent if one 
differentiates Eq. 9 with respect to time, giving a sec- 
ond order, harmonic oscillator-like equation, plus a 
dissipativc coupling to the bulk concentration. 

Whether the surface concentration or the bulk con- 
centration near the surface vanishes at long times de- 
pends on the surface kinetics as well as geometry. If 
the combined surface kinetics towards annihilation is 
slow relative to bulk diffusion, the decay of the sub- 
layer concentration n s can be extremely slow. Simi- 
larly, if the number of surface states is large, there is 
an effective delay to annihilation and a higher proba- 
bility that a surface species can detach and replenish 
the sublayer concentration. This effect is very sen- 
sitive to the annihilation rate p* and the size of the 
reaction N and can keep the subsurface concentration 
high for essentially all times. 

A number of extensions and related approaches 
to this and related systems can be readily investi- 
gated. For example, in applications such as sur- 
factant adsorption, the surface concentration T can 
be appreciable and suppress additional adsorption. 
If surface species Ti has molecular area eij, an ad- 
sorption term including steric exclusion would be 
jn s (l — X^iLi ^;)- The surface rate equations re- 
main linear in T, but with a time-dependent tran- 
sition matrix M. The effective boundary condition 
Eq. 9 is now nonlinear in n s (t) through -F[n s ]. How- 
ever, for many biochemical applications (such as cell 
signaling and virus adsorption and entry) the total 
surface concentration is low such that X) i=i a i^i 1 
and the adsorption term can be linearized. In our 
one-dimensional analysis, as long as N is not too large 
and there is an appreciable annihilation process, the 
surface concentrations all vanish in time. 

The effects of multistage adsorption can also be 
explored on surfaces of arbitrary shape, particularly 
for cylinders and spheres. For multistage processes 
on a sphere, the sublayer concentration approaches a 
positive value n s (t — > oo) = n (l — l ^° Kq )■ We also 



expect positive eigenvalues Xj > of M to have a 
striking effect on the transport. 

Finally, although we have only considered sim- 
ple linearized surface reaction schemes with negative 
eigenvalues, systems that support oscillations, such 
as those involved in surface-mediated cell signalling, 
could also be treated within our framework. Features 
of the surface reactions and the bulk concentrations 
near the reacting surface remain coupled through 
the kernel K(t). Under certain conditions, nonlin- 
ear surface reaction schemes may also be linearized. 
One example is in the stochastic representation of 
the surface reactions. If we write the surface quanti- 
ties in terms of the probability distribution function 
P(ni, 7i2, ^3, • • • , t) that there are n\ molecules of of 
type 1, ri2 of type 2, etc., the surface reactions can 
be written as a linear Master equation. This allows 
our approach to be applied when Ti in the last term 
of Eq. 4 is interpreted as (ri(£)), the ensemble av- 
erage X){n } n iP({ n i}i t). Using this interpretation, 
the full problem can be solved using linear methods 
similar to those presented, albeit for extremely large 
matrix dimension. 
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VI. APPENDIX 

Here, we derive asymptotic expression for bulk and 
surface densities. The trivial short time behavior of 
the subsurface density is n s (t)/no ~ (1—jt), indepen- 
dent of the surface reactions since the first physical 
phenomenon to occur is particle adsorption from the 
bulk to the interface, at rate 7. 

For large distances (z — l)/y/~Dt 3> 1 and in the 
limit jK(s = 0) = jKq y/D/t, asymptotic eval- 
uation of the inverse Laplace inversion integral over 
n(z, s) yields 



n(z,t) 
n 



1- 1 



lK Vt 



exp 



{z-lf 



ADt 

"(20) 

The condition ^/Kq 3> y/D/t can be interpreted as a 
comparison between two typical velocities. The usual 
diffusive velocity, \J D/t, is compared to an effective 
reaction velocity expressed by 7 modulated by surface 
effects through the kernel Kq. We may thus define 
an effective Damkohlcr number D a = jKoy/t/^/D, so 
that Eqn. 20 is valid only at large distances and for 
large values of D a . The leading term on the right- 
hand-side above is independent of the surface kinet- 
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ics: the first information to have traveled away from 
the interface is the initial depletion of the n s layer 
into the surface and interfacial effects emerge as first 
order corrections. 

In the t — > oo limit, the dominant contribution to 
n(z,t) comes from small values of u in Eq. 16. Ap- 
proximating L(l,u) with its u — > limit, we find the 
asymptotic long time limit 




no 



D 



(21) 

This expression is valid only if the surface dynamics 



include a net sink of material. As long as there is 
some annihilation, K = 1 + Y$Li V^Vji/Xj < 
and Eq. 21 holds. A similar consideration of the 
small- u dominated integration in Eq. 17 yields for the 
surface concentrations 



r*(|V|( 



where A* < is the largest eigenvalue of the chemi- 
cal transition matrix M. Both Eq. 21 and 22 exhibit 
diffusion-limited \j\ft behavior. These general re- 
sults rely only on the linearity of -F[n s ] and are valid 
for any surface reaction scheme through the eigenval- 
ues and eigenvectors of the transition matrix M and 
the resulting function K(s). 
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